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ABSTRACT 

N-body + hydrodynamic simulations of galaxy clusters are used to demonstrate a correlation between 
galaxy cluster mass and the strength of the Sunyaev-Zel'dovich (SZ) effect induced by the cluster. 
The intrinsic scatter in the correlaton is larger than seen in the cluster mass — X-ray temperature 
correlation, but smaller than seen in the correlation between mass (or temperature) and X-ray luminosity, 
as expected. Using the convergence to self-similarity of cluster structure at larger radii, a simple area- 
averaged SZ value derived from mock SZ maps also correlates well with mass; the intrinsic scatter in 
this correlation is comparable to that seen in simulations for the mass — temperature correlation. Such 
a relation may prove a powerful tool for estimating cluster masses at higher redshifts. 

Subject headings: Galaxies-clusters, cosmology-theory 
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1. MOTIVATION 

Correlations between galaxy cluster mass and observ- 
ables such as X-ray luminosity L x or temperature T are 
used extensively to probe cosmology. For instance, the 
abundance of clusters as a function of temperature has 
been used to constrain the normalization of the power 
spectrum of primordial fluctuations, a^, while searches for 
evolution in the cluster X-ray luminosity function have 
been employed to argue for values of the cosmological den- 
sity param eter fin. Theoretica l models such as the Press- 
Schcchter ( Press k Schechtei 1974) formalism character- 
istically provide the abundance of clusters as a function 
of the cluster mass. The use of temperature or luminos- 
ity functions to constrain the parameters in the theoret- 
ical model therefore depends upon a relation between T 
or L x and mass. Theoretical arguments and cluster sim- 
ulations indicate that temperature and X-ray luminosity 
should have a power law dependence on mass, while ob- 
servations of cluster abundances and correlations between 
observables are consistent with such simple power law re- 
lations holding true for galaxy clusters. 

A comparatively new tool for investigating clusters is 
provid ed by observations of t he Su nyaev-Zel'do vich (SZ) 
effect flSunyaev k Zel'dovichl 1972; |Birkinshaw| 1998) - 



the distortion in the cosmic microwave background spec- 
trum produced by Compton upscattering of microwave 
background photons by electrons in the hot intraclus- 
ter medium (ICM). The SZ effect has two properties 
that make it particularly interesting for examining high- 
rcdshift systems. First, the signal strength is not attenu- 
ated by redshift, in contrast with X-ray flux. Second, for 
clusters of fixed mass, self-similar scaling relations predict 
a strong positive evolution with redshift. Since observa- 
tions of X-ray temperature, X-ray surface brightness, and 
the SZ effect depend upon the density and temperature 
structure of cluster gas in different ways, the SZ effect 
gives a complimentary probe of the state of the ICM. 

An interesting question is whether the strength of the 



SZ effect correlates with cluster mass, as is expected of 
T or L x . A simple argument reproduced below suggests 
that such a correlation should also be present; verifica- 
tion of this would provide a consistency check upon the 
simple theoretical models used to describe clusters. Fur- 
thermore, accurate X-ray temperature determinations be- 
come harder with increasing redshift. As the strength of 
the SZ signal is redshift-independent, and as clusters of a 
fixed mass produce a stronger effect at higher redshift, a 
correlation between SZ signal and mass yields a powerful 
method for estimating cluster masses at redshifts beyond 
those currently probed by X-ray temperature estimates. 

In this Letter, we use simulated clusters to examine the 
relation between cluster mass and the SZ effect. We con- 
struct mock SZ images of the simulated clusters, extract 
SZ observables, and examine correlations with mass as 
well as scatter about the correlations. We find that the 
central SZ signal exhibits the expected self-similar scaling 
with mass in these simulations, with a scatter significantly 
smaller than is seen in simulated L x -M relations or the 
observed cluster L x -T relation. We will also note that a 
simple area-averaging of the SZ signal reduces the disper- 
sion to a level near that expected from the cluster X-ray 
temperature. Some concerns about applying this result, 
as well as directions for future study, will be described at 
the end. 

2. MODELS 

The simulated clusters were produced using the N-body 
+ hydrodynamic code P3MSPH; both the simulati on code 
and analysis techniques are described in Evrarcj (1990). 
Each simulation run produced an individual cluster, with 
initial conditions drawn from the f2 = 1, f2& = 0.1, 
h = 0.5 (assumed throughout this paper), as = 0.6 
cold dark matter model using the constrained-realization 
technique of Bcrtschinger (1987). The clusters simulated 
ranged in mass from 8 x 10 13 to 1.1 x 10 15 M Q , with ICM 
temperatures between 0.7 and 4.3 keV. A total of 65536 
particles were used in each simulation, half in dark matter 
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and half in baryons. The fractional mass resolution was 
comparable amongst the simulated clusters, with around 
3500 gas particles typically residing in the cluster at z = 0. 
A total of 73 simulated clusters were produced in this 
fashion. The final, z = configuration of each simu- 
lated cluster was then "imaged" in the SZ by producing 
a pixel map of values of the Compton y-parameter, each 
pixel's value constructed by a line-of-sight pressure inte- 
gral through the simulated cluster. For each run, an SZ 
map was produced in three orthogonal directions; this pro- 
duced a total of 219 maps. In extracting a central value 
for the y-parameter, y , the cluster "center" is defined by 
the projected location of the bottom of the cluster poten- 
tial well, which traces well the projected location of the 
maximum value of the y-parameter. 

3. THE CLUSTER Y-M RELATION 

The Compton y-parameter induced by upscattering of 
microwave background photons along a line of sight I 
through the cluster is given by 
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(1) 



This equation can be used to derive a self-similar scal- 
ing relation for clusters which are near isothermality when 
density- weighted. Writing the density in terms of the 
background density p, expressing length scales in terms 
of the radius r$ at some fixed overdensity 5, and factoring 
out the temperature, we have 



y oc p r s T 
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If clusters are self-similar, then the integral produces 
a number which is ident ical f or clust ers of any mass 
dNavarro, Frenk & White] 1995; |Metzlerj 1995; [Metzler & 
Evrard 1998), and varies only with the choice of line of 



sight. With M = Sp (4?r/3) rf, p oc (1 + z) , and the 
virial scaling T oc M 2 / 3 (1 + z), we have 



y (x M(l + zy 



(3) 



This yields the scaling, with mass and redshift, of the y- 
parameter induced along a line of sight with a 2-D radius 
equal to the 3-D radius of some fixed fiducial overdensity; 
the constant of proportionality depends upon the overden- 
sity in question. In particular, the values of the SZ central 
temperature decrement should scale in this fashion. 

Do the simulated clusters obey this relation? We have 
219 pairs of yo and mass within an overdensity of 200 with 
which to test the relation; the results are plotted in Figure 
1. The best-fit power law relation is plotted as a solid line; 
the best-fit power law slope is 0.98 ± 0.04, well in agree- 
ment with the theoretical expectation, Eq. (||). Demand- 
ing Eq. (||) be exactly correct yields a best-fit coefficient 
of proportionality which is not significantly different from 
that found in the general power-law fit. 

While the success of the expected scaling law is quite 
striking, the small scatter about the mean relation is also 
of note. Modelling the distribution of possible values of the 




(M,) 



FlG. 1. — Mass within an overdensity of 200 versus central y- 
distortion, for the 219 simulated cluster observations in the datasct. 
Each point marks an observation. A solid line marks the least-squares 

best-fit to a power law, y = lO" 4 - 4 * - 7 (Af/1O 15 M )°' 98±0 ' 04 . A 
linear relation, expected by theory, is not significantly different; forcing 
a fit to a linear relation produces the same coefficient, and the dashed 
line on the figure marking this relation cannot be distinguished from the 
solid line. 



y-parameter for a cluster of a given mass as log-normal, 
the simulated data show a dispersion of a = 0.15 in log yo- 
While larger than the scatter predicted by simulations 
for the T -M re lation (Evrard, Metzler & Navarro 1996, 
hereafter EMN), this dispersion is significantly smaller 
than the scatter predicted b y the lu minosity-mass rela- 
tion observed in simulations ( Mctzlerj 1995), and observed 
in th e luminosity-tempe rature correlation for X-ray clus- 
ters ( |Arnaud fc Evrard| 1998). 



The magnitude of the intrinsic scatter about the mean 
SZ-mass relation is important. Since the predicted abun- 
dance of clusters falls with increasing mass, the net effect 
of such scatter is to artificially populate the "bright" end 
of the abundance function. In other words, such an intrin- 
sic scatter would make the real universe appear to have 
more high-mass clusters than are actually present. Since 
the high-mass end of the the cluster abundance function 
also contains most of the cosmological constraining power, 
a large intrinsic dispersion in the correlation of an observ- 
able with mass makes using the cluster abundance to con- 
strain cosmological parameters extremely difficult. Most 
efforts to construct the cluster abundance include an at- 
tempt to correct for this. However, an accurate correction 
requires an understanding of the intrinsic dispersion in the 
correlation with mass, which has not been available. 

Consequently, the comparatively small scatter suggests 
that the SZ central decrement is a much more robust ob- 
servable to use in constructing an abundance function than 
X-ray luminosity. If this result holds at higher redshift as 
well, then along with the redshift-independence of the SZ 
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signal, and the expected self-similar scaling of the SZ sig- 
nal strength with redshift at fixed mass, we can expect a 
large dataset of measurements of the SZ central temper- 
ature decrement to yield a much better estimate of the 
high-redshift cluster abundance. 

4. THE CLUSTER (Y) A -M RELATION 

Cluster entropy profiles drawn from simulations such 
as used here typically exhibit self-similar behavior out- 
side of the cluster core, but a dispersion in scaled en- 
tropy values in the center of the cluster. Observations 
suggest such a dispersion in central entropy states, as we 
can demonstrate. First, note that the cluster luminosity- 
temperature relation is both theoretically and empirically 
well-described by a power law, L x oc T Q . Second, 
the fact that most emission originates from the cluster 
core allows us to approximate the X-ray luminosity by 
L x oc n 2 , T 1 / 2 rj!, where n is the central gas density and 
r c is the cluster X-ray core radius. Equating, we can solve 
for the quantity s' = T/rt^ 3 , which is simply related to 
the entropy by s = logs'; we obtain 



s cx r. 



rp(7-2a)/6 



(4) 



If the relation between central entropy and cluster mass 
had little intrinsic scatter, then the tight mass — temperature 
relation would imply a small dispersion in observed values 
of cluster X-ray core radius at a given X-ray tempera- 
ture. In other words, we should see a tight relation be- 
tween X-ray core radius and X-ra y tem perature; this is 
contradicted by observations ( Mohi 1998; [Mushotzky 1998 



1998). At fixed temperature, thus mass, there is a wide 
range of observed core radii, and thus by inference central 
entropy states. 

Therefore, the result that the magnitude of the scatter 
in the cluster y — M relation lies between that of the T- 
M and L x -M relations is not s urprisi ng. While the T-M 
relation is expected to be tight (EMN), as the gas specific 
energy is determined by the potential the gas falls through, 
the dispersion in central entropy noted above results in a 
dispersion in central gas density values. The y-parameter 
depends upon an integral of the gas density along the line 
of sight, so such an intrinsic dispersion introduces scatter 
in the correlation with mass. However, the X-ray emission 
depends upon an integral of the square of the density; 
the intrinsic scatter in the luminosity — mass relation is 
correspondingly larger. 

That the gas density (and thus gas entropy) appears in 
simulations to converge to self-similar behavior at larger 
radii suggests that the scatter in correlations between mass 
and ICM-related properties could be reduced through the 
use of a statistic which increases the weight of data from 
outside the cluster core. For X-ray emission, this has been 
demonstrated through the prediction and observation of 
a tight correlation between cluster temperature a nd ra- 
dius conta ining a fixed amount of X-ray emission (Mohi 
fc Evrard 1997). It should be possible to construct an 
SZ observable which correlates well with mass and shows 
much smaller intrinsic disperson at fixed mass than the 
central y-distortion. 



Defining our new statistic as {u)a, we can write {h)a 
in terms of the y-distortion induced by a line-of-sight at 
angle 9 from the cluster center by 



y(6)W(9)d 2 



(5) 



where W (6) is a window function that determines the rel- 
ative weighting of different parts of the cluster. The sim- 
plest approach is to use a top-hat window function; {y) a is 
then simply an average of the SZ map within some chosen 
radius. 

It is important that the radius used to construct {y)A 
be a self-similar radius, i.e. the radius out to some fidu- 
cial overdensity S in 3-D, r$. This ensures that the same 
regions of different clusters are being examined in con- 
structing (ij)a- In practice, however, we would not know 
the radius r$ for an observed cluster. It could be deduced 
from the X-ray temperature; but if we have a good mea- 
surement of the temperature, then we immediately have 
access to a mass estimator which is independent of the gas 
density, and thus presumably has smaller intrinsic error. 
The point is to develop an SZ-based mass estimator to 
be used at higher redshifts, where X-ray temperatures are 
difficult to obtain. 

We solve this problem by estimating r$ using the cen- 
tral value of the y-parameter alone. In the previous 
section, we demonstrated a correlation between the cen- 
tral y-distortion and cluster mass at r2oo! we use this 
relation to estimate r200j then use this value to set an 
outer limit for averaging an SZ map to construct {y)A- 
Here we will use 0.3r2oo, corresponding to 5 ~ 2000 
for the NFW density profile. With our previous result 
of y = 3.98 x lO^Af^ 98 , where M 15 is the mass 
within an overdensity of 200 in units of 1O 15 M0, and 
with Monti 



200 p x 47rr| 00 /3 by definition, we obtain 



1.61 Q, 



-1/3 



yo/io- 



.1/3 



Mpc. For each mock im- 



age, we extract yo, use this value to find t^qq, and then con- 
struct the average {y)A within 0.3r2oo; this radius would 
correspond to 6 « 1.26' for yo = 10~ 5 at z = 0.5. 

The result is shown in Figure 2, with a best-fit power 
law shown as a solid line, while the best fit to Eq. (^) is 
shown as a dashed line. The two are barely distinguish- 
able. The slope of the best fit relation is 0.97±0.01, slightly 
shallower than the expected slope; the plot shows, how- 
ever, that a linear scaling provides a good fit as well. More 
notably, the dispersion in the relation is much reduced; 
when modelled as log-normal, the scatter in log(y)A at 
fixed mass is only 0.06 — that is, 13 — 15% in (y)A — much 
smaller than for the central decrement alone and compa- 
rable to that found by EMN for the mass-temperature 
relation at overdensities of ^500. 

It may seem puzzling that the scatter in the \og{y) a ~ 
M relation is less than that in the yo — M relation, given 
that we use yo to choose a radius within which to construct 
\og(y)A- However, note first that the dispersion in r|g will 
be only 1/3 the dispersion in mass inferred from the cen- 
tral signal. Second, the fact that the SZ signal should fall 
with cluster radius suppresses the error introduced by us- 
ing an incorrect inferred r2oo- For instance, imagine that 
a cluster with a given mass has an observed yo higher than 
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Fig. 2.- 



Mass within an overdensity of 200, M20O1 ver_ 



sus average y-distortion over a circular area of radius 0.3 r^Q, 



(y)A- 



r|jjQ is an estimate of the three-dimensional radius con- 
taining an overdensity of 200, derived from the central SZ signal 
yo. Each point marks one of the 219 simulated cluster observa- 
tions. A solid line marks the least-squares best-fit to a power law, 

- , , „ „ / , , - \ 0.97±0.01 
(y) A = iq-5-1±0-3 (M/10 15 Mq) . Forcing a fit to a lin- 

ear relation produces the same coefficient; this curve is drawn in a 
dashed line on the figure. 



expected from the mean relation; t^qq will then be higher 
than the correct value. However, since the signal falls with 
radius, constructing an average within this incorrect radius 
depresses the result, counteracting the impact of the over- 
large central value. This works to suppress the dispersion 
in (y)A for clusters at a fixed mass that would be induced 
by th e dispersion in yg. 



5. DISCUSSION 

The simulation results discussed above evidence a cor- 
relation between the strength of SZ signal and collapsed 



in agreement with theoretical expectation. The in - 
trinsi d scatter observed in the relation between the central 
value pi the Compton }r-paramcter yp and mass is greater 



than aeen in the mass — temperature relation , hnt consiHer- 
ably smaller than simulations show in the relation between 
X-ray luminosity and mass, as expected. The scatter in 
this relation can be reduced to a level near that of the mass 
— temperature relation through the use of statistics which 
depend less sensitively on central gas and more sensitively 
on gas at larger radii. 

Ssomc cautionary notes are in order. First of all, the op- 
timal statistic to use will in all likelihood be experiment- 
dependent, since different experiments sample the sky dif- 
ferently. Also, such a correlation with mass should be af- 
fected by limits in observational resolution; when a cluster 
is not resolved, the y-parameter measured will scale with 
the angular-diameter distance to the cluster as y cx d^ 2 . 



There are also issues of concern about the simulations 
used here. Chief among these are numerical resolution lim- 
itations, which likely serve to reduce the dispersion in the 
SZ — mass correlation. Furthermore, the effect of miss- 
ing physics which breaks self-similarity, such as cooling or 
the inclusion of supernova-driven galactic winds, must be 
considered. Next, the simulated clusters used in this inves- 
tigation were taken at z = 0; numerical resolution issues 
are of greater concern as redshift increases. It would be 
worthwhile to investigate in detail, using higher-resolution 
simulations, the status of the SZ — mass correlation at 
redshifts of 0.5 or 1. Finally, the effects of the choice of 
background cosmology — particularly, the values chosen 
for f^o and the baryon fraction = f^/Slo — must be 
examined. 

One hope for verifying or calibrating the SZ — mass 
correlation comes from comparing SZ-derived masses with 
those from X-ray observations or from gravitational lens- 
ing. Since low redshift clusters typically have better tem- 
perature estimates, the planned Viper Sunyaev-Zel'dovich 
Survey (FLomer 1998) should provide a good opportunity 
to test this correlation. 

Future work is planned to study these issues, as well 
as to examine the systematics of using SZ observations to 
constrain the cluster baryon fraction. 
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